5 Path Integrals in Quantum Mechanics and Quan- 
tum Field Theory 

In the past chapter we gave a summary of the Hilbert space picture of Quantum 
Mechanics and of a Scalar Quantum Field Theory. Here we will present the Path 
Integral picture of Quantum Mechanics ans Scalar Quantum Field Theory. 

The Path Integral picture is important for two reasons. First, it offers an 
alternative, complementary, picture of Quantum Mechanics in which the role 
of the classical limit is apparent. Secondly, it gives a direct route to the study 
regimes where perturbation theory is either inadequate or fails completely. A 
standard approach to these problems is the WKB approximation (Wcntzel, 
Kramers and Brillouin). As it happens, it is extremely difficult (if not im- 
possible) the generalize the WKB approximation to a Quantum Field Theory. 
Instead, the non-perturbative treatment of the Feynman path integral, which 
is equivalent to WKB, is generalizable to non-perturbative problems in QFT. 
In this chapter we will use path integrals only for bosonic systems, such as 
scalar and abelian gauge fields. In the following chapters we will also give a 
full treatment of the path integral, including its applications to fermionic fields, 
non-relativistic many body systems, and Abelian gauge fields. Non-Abelian 
gauge fields will be discussed in the following course, Physics 583. 

5.1 Path Integrals and Quantum Mechanics 

Consider a simple quantum mechanical system whose dynamics can be described 
by a (generalized) time-dependent coordinate operator q(t), i.e., the position 
operator in the Heisenberg representation. We will denote by \q, t) an eigenstate 
of q{t) with eigenvalue q, 

q(t)\q,t) = q\q,t) (1) 
We want to compute the amplitude 

F{q f ,t f \qi,ti) = (q f ,t f \q h ti) (2) 

Let qs be the Schrodinger operator, related to the Heisenberg operator q(t) by 
the action of the time evolution operator: 

q(t) = e lAt l h q s (3) 
and whose eigenstates are \q), 

qs\q) = q\q) (4) 
The states \q) and \q,t) are related via the evolution operator 

\q)=e- i ^ h \q,t) (5) 

Therefore the amplitude F(qf ,tf\qi,ti) is a matrix element of the evolution 
operator 

F(q f ,t f \ qi ,t t )^(q f \e^-^^\q t ) (6) 
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The amplitude F(qf,tf\qi,ti) has a simple physical interpretation. Let us set, 
for simplicity, \qi,ti) = |0, 0) and \q/,tf) = \q,t). Then, from the definition of 
this matrix element, we find out that it obeys 

|imF(g,i|0,0) = <#>=*(<?) (7) 
Furthermore, after some algebra we also find that 
OF d 

= m^- t (q\e- i6t / h \0) 

= {q\He- ikt ' h \Q) 

= J dq'iqlHWjiq'le-^O) 



(8) 



where we have used that, since {\q}} is a complete set of states, the identity 
operator / has the expansion, called the resolution of the identity, 

1 = J dq'\q'){q'\ (9) 

Here we have assumed that the states are orthonormal, i.e., 

(q\q')=S(q-q') (10) 

Hence, 

dF f 

ih— = I dq'(q\H\q')F(q', t\0, 0) = H q F(q, t\0, 0) (11) 

In other terms, F(q 7 t\0,0) is the solution of the Schrodinger Equation that 
satisfies the initial condition of Eq. 7. The amplitude F(q, t\0, 0) is called the 
Schrodinger Propagator. Let us next define a partition of the time interval 
[U,tf] into N intervals of length At, 

t f -ti = NAt (12) 

Let {tj} (with j = 0, .. ., N + 1) denote a set of points in the interval [t,tf], 
such that 

U — to < ti < . . . < ijy < tN+i — tf (13) 

Clearly, t k = t + kAt (for k = 1, . . . , N + 1. 
The superposition principle tells us that 

F(q f ,t f \q i ,t i ) = J dq'{q f ,t f \q',t l ){q',t'\q i ,t i ) (14) 
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where we have inserted I in the form of the resolution of the identity of Eq. 14. 
By repeating this process many times we find 

F(qf,tf\qi,U) = Jdq 1 ...dq N (qf,t f \q N ,t N )(q N ,t N \q N ^ 1 ,t N ^ 1 ) x... 

x . . . (q j ,t j \q j - 1 ,t j - 1 ) . . . (q 1 ,t 1 \q l ,t i ) 

(15) 

Each factor (qj, tj\qj-i, tj-i) has the form 

(^.t^-l^-l) = (q,\e-^^/%^) EE {q^^q^) (16) 

In the limit N — ► oo, with \tf — U\ fixed and finite, the interval At becomes 
infinitesimally small and At — > 0. Hence, as N — > oo we can write the approxi- 
mate expression for (qj, tj\qj-\, tj-i), 

(QiMU-uti-i) = fele-^/^-i) 

= (qjlI-i^H + OdAtf)^) 



= 6( qj qj.,) - l ^-(q 3 \H\q 3 _ 1 ) + Q{(Atf ) 



(17) 



which becomes asymptotically exact as ./V — > oo. 

We can also introduce a complete set of momentum eigenstates {\p}} and 
their resolution of the identity 

/oo 
dp\p)(p\ (18) 

Recall that the overlap between \q) and \p) is 

(fl|p> = -2= eta/ft (19) 
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Figure 1: A history q(t) of the system. 



For a typical Hamiltonian of the form 



its matrix elements are 



)/h 



(20) 



(21) 



Within the same level of approximation we can also write 



(Qj,tj\lj-i,tj-i) ~ y ^" ex P 



(22) 

where we have introduce the "mid-point rule" which amounts to the replacement 
1j ~ * \ (ij + inside the Hamiltonian H(p, q). Putting everything together 

we find that the matrix element (qf,tf\qi,ti) becomes 



r N „oo N+l 

(q f ,t f \qi,ti) = lim JJdqj JJ 

J —CO j — l 



. N+l r 



27rft 



Pjilj - - At H [Pj 



1i + 



(23) 



Therefore, in the limit N — > oo, holding |i — 1/| fixed, the amplitude (qf 7 tf\qi, ti) 
is given by the (formal) expression 



(qf,t f \qi,ti) = J VpVq eJ 



dt [pq - H(p, q)} 



(24) 
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where we have used the notation 



VpVq= lim TX^A (25) 

i=i 

which defines the integration measure. The functions, or configurations, (q(t),p(t)) 
must satisfy the initial and final conditions 

q(U)=qi, q(t f )=q f (26) 

Thus the matrix clement (qf,tf\qi,ti) is expressible as a sum over histories in 
phase space. The weight of each history is the exponential factor of Eq. 24. 
Notice that the quantity in brackets it is just the Lagrangian 

L=pq-H(p,q) (27) 

Thus the matrix element is just 

(q f ,t f \q,t)= J VpVqe iS ^' h (28) 

where S(q,p) is the action of each history (q(t),p(t)). Also notice that the sum 
(or integral) runs over independent functions q(t) andp(t) which are not required 
to satisfy any constraint (apart from the initial and final conditions) and, in 
particular they are not the solution of the equations of motion. Expressions of 
these type are known as path-integrals. They are also called functional integrals, 
since the integration measure is a sum is over functions, instead of numbers as 
in a conventional integral. 

Using a Gaussian integral of the form (which involves an analytic continua- 
tion) 

p 2 . At , At , 

f^ e ^-^T=./^Z e ^ 2 (29) 

we can integrate out explicitly the momenta in the path-integral and find a 
formula that involves only the histories of the coordinate alone. (Notice that 
there arc no initial and final conditions on the momenta since the initial and 
final states have well defined positions.) The result is 

r z / dt L(q, q) 
(q f ,t f \q u U) = JVq e h Ju (30) 

which is known as the Feynman Path Integral. Here L(q, q) is 

L(q,q) = ±mq 2 -V(q) (31) 

and the sum over histories q(t) is restricted by the boundary conditions q(U) = qi 
and q(t f ) = q f . 



5 



q(tf) 



Figure 2: Two histories with the same initial and final states. 



Notice now that the Feynman path-integral tells us that in the correspon- 
dence limit H — > 0, the action S must be stationary since otherwise the con- 
tributions of the rapidly oscillating exponential will add up to zero. In other 
words, in the classical limit there is only one history q c (t) that contributes. For 
this history, q c (t), the action S is stationary, SS = 0, and q c (t) is the solution 
of the Classical Equation of Motion 

Q ± **± (32) 

In other terms, in the correspondence limit K — > 0, the evaluation of the Feynman 
path integral reduces to the requirement that the Least Action Principle should 
hold. We are in the classical limit. 

5.2 Evaluating Path Integrals in Quantum Mechanics 

Let us first discuss the following problem. We wish to know how to compute 
the amplitude (qf,tf\qi,U) for a dynamical system whose Lagrangian has the 
standard form of Eq. 31. For simplicity we will begin with a linear harmonic 
oscillator. 

The Hamiltonian for a linear harmonic oscillator is 

H= k + —* (33) 

and the associated Lagrangian is 



m 



L = jq 2 - —q 2 (34) 



G 



Let q c (t) be the classical trajectory. It is the solution of the classical equations 
of motion 



d 2 ' 



dt 2 



+ u 2 q c 



= 



(35) 



Let us denote by q(t) an arbitrary history of the system and by £(t) its deviation 
from the classical solution q c (t). Since all the histories, including the classical 
trajectory q c (t), obey the same initial and final conditions 

q(U) = % q(tf) = Qf (36) 

it follows that obeys instead vanishing initial and final conditions: 

ttU) = £(*/) = (37) 

After some trivial algebra it is easy to show that the action S for an arbitrary 
history q(t) has the property 



S{q,q) = S{q c ,q c ) + S{Z,i) + 



dt 



dt 



dq c 
dt 



+ 



dt 77l£ 



'<lc 



dt 2 



+ uj 2 q c 



(38) 

The third term vanishes due to the boundary conditions on £, Eq. 37, and the 
last term also vanishes since q c is a solution of the classical equation of motion 
Eq. 35. These two features hold for all systems, even if they are not harmonic. 
However, the Lagrangian (and hence the action) for £, the second term in Eq. 
38, in general is not the same as the action for the classical trajectory (the first 
term). Only for the harmonic oscillator S(t;, £) has the same form as S(q c ,q c ). 
Hence, for a harmonic oscillator we get 



(qf,t f \qi,U) = e 



J 



€(*/)=0 



(39) 



Notice that the information on the initial and final states enters only through 
the factor associated with the classical trajectory. For the linear harmonic 
oscillator, the quantum mechanical contribution is independent of the initial 
and final states. Thus, we need to do two things: 1) we need an explicit solution 
q c (t) of the equation of motion, for which we will compute S(q c ,q c ), and 2) we 
need to compute the quantum mechanical correction, the last factor in Eq. 39, 
which measures the strength of the quantum fluctuations. 

For a general dynamical system, whose Lagrangian has the form of Eq. 31, 
the action of Eq. 38 takes the form 



S(q,q) = S(q c ,q c ) + S e ti(€,£;q c ) 



/ d 2 q c 8V 



(40) 
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where the effective Lagrangian L e ff(£, £; <Zc) is 

c 2 - o(e) (4i) 



ieff(e,C) = ^m^ 2 -- 



2 s 2 a? 2 



Qc 



Once again, the boundary conditions £(U) = £(t/) = and the fact the q c (t) is 
a solution of the equation of motion together imply that the last two terms of 
Eq. 40 vanish identically. 

Thus, to the extent that we are allowed to neglect 0(£ 3 ) the corrections, 
the effective Lagrangian i off can be approximated by a Lagrangian which is 
quadratic in the fluctuation £. In general, the effective Lagrangian will depend 
on the actual classical trajectory, since V (q c ) in general is not a constant, 
but a function of time determined by q c (t). However, if one is interested in 
the quantum fluctuations about a minimum of the potential V(q), then q c {t) is 
constant (and equal to the minimum). We will discuss below this case in detail. 

Before we embark in an actual computation it is worthwhile to ask when 
should the neglecting of the 0(£ 3 ) terms be a good approximations. Since we 
are expanding about the classical path q c , we expect that this approximation 
should be correct as we formally take the limit h — > 0. In the path integral 
the effective action always appears in the combination S c s/h. Hence, for an 
effective action which is quadratic in £, we can eliminate the dependence on h 
by the rescaling 

£ = (42) 
This rescaling leaves the classical contribution S(q c )/h unaffected. However, 
terms with higher powers in £, say 0(£"), scale like h n ^ 2 . Thus the action 
(divided by K) has an expansion of the form 

| = ^(q c ) + S^(i;q c ) + J2^ /2 S (n) (l,qc) (43) 

n— 3 

Thus, in the limit h — ► 0, we can expand the weight of the path integral in 
powers of h. The matrix clement we are calculating then takes the form 

(q^tf^^^e^^^ Z^{q c ) (l + 0(h)) (44) 

The quantity (q c ) is the result of keeping only the quadratic approximation. 
The higher order terms are a power series expansion in h and are analytic in h. 
(Here I have used the fact that by symmetry the odd powers in £ in general do 
not contribute, although there are some cases where they do.) 

Let us now calculate the effect of the quantum fluctuations to quadratic 
order. This is the WKB approximation. Let us denote this factor by Z, 

Z (2) {qc) = [ ^ J Sl! dt L&> (llqc) (45) 

It is elementary to show that, due to the boundary conditions, the action 
£eff(£>£) becomes 



sMi.i) = \J ' di id) 



■2 



- W ^2 -V "(&(*)) 



m (46) 
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The differential operator 



A=-m^-V"(q c (t)) (47) 



has the form of a Schrodinger operator for a particle on a "coordinate" t in a 
potential — V (q c (t)). 

Let il> n (t) be a complete set of eigenfunctions of A satisfying the boundary 
conditions ip(ti) — ip(tf) = 0. Completeness and othonormality implies that the 
eigenfunctions {ip n (t)} satisfy 

n 

/ dt ip* n {t)\jj m (t) = 5 nim 

(48) 

We can expand £(t) as a linear combination of the eigenfunctions {tp n (t)}, 

n 

Clearly, we have = — as we should. 

For the special case of qi — qf — qo, where go is a minimum of the potential 
V(q), V"(qo) = cjcff > is a constant, and the eigenvectors of the Schrodinger 
operator are just plane waves. (For a linear harmonic oscillator c<j c ff = uj.) Thus, 
in this case the eigenvectors are 

Mt) = b n sm(k n (t-U)) (50) 

where 

7TTI 

k n = - n= 1,2,3,... (51) 

t f - U 

and b n = 1/ ' y/tf — U. The eigenvalues of A are 

A n = k 2 n - u; 2 s = j——^n 2 - u; 2 s (52) 
By using the expansion of Eq. 49, we find that the action takes the form 



(53) 



where we have used the completeness and orthonormality of the basis functions 
{tpn(t)}- Furthermore, the expansion is a canonical transformation £(i) — > c„. 
More to the point, the expansion is actually a parametrization of the possible 
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histories in terms of a set of orthonormal functions, and it can be used to define 
the integration measure to be 



n v 

with unit Jacobian. Here Af is an irrelevant normalization constant that will be 
defined below. 

Finally, the (formal) Gaussian integral (defined by a suitable analytic con- 
tinuation) 



r 

j — c 



d% e (i/2)A n cl = [_ iAi ]-V2 (55) 



/27T 

can be used to write the amplitude as 

=N JjA- 1 / 2 = Af (DetA)- 1 / 2 (56) 

n 

where we have used the definition that the determinant of an operator is equal 
to the product of its eigenvalues. Therefore, up to a normalization constant, 
Z^ = (DetA) -1 / 2 . We have thus reduced the problem of the computation of 
the leading (Gaussian) fluctuations to the path-integral to the computation of 
a determinant of the fluctuation operator, a differential operator defined by the 
choice of classical trajectory. Below you will see how this is done. 

It is useful to consider the related problem obtained by an analytic continua- 
tion to imaginary time, t — ► —it. We saw before that there is a relation between 
this situation and Statistical Physics. We will now work out one example that 
will be very instructive. 

Formally, upon the analytic continuation t — ► —it we get 

(q f \ exp ^H(t f - ti)j \ qi ) - (q f \ exp [~\h{t } - r^j (57) 
Let us choose 

n = T f =/3h (58) 

where = 1/T, and T is the temperature (in units of ks = 1)- Hence, we find 
that 

(q f ,-i0/h\ qi ,O) - (q f \e-0 H \ qi ) (59) 

The operator p 

p = e~ 0H (60) 

is the Density Matrix in the Canonical Ensemble of Statistical Mechanics for a 
system with Hamiltonian H in thermal equilibrium at temperature T. 
It is customary to define the Partition Function Z, 



Z = tie-' m = J dq (q\e- 0H \q) 



(61) 
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where I inserted a complete set of eigenstates of Q. Using the results that were 
derived above, we see that the partition function Z can be written as a Feynman 
path integral in imaginary time, of the form 




0h 



dr 



1 (dq 
—m — — 

2 \dr 



-Jo"' ^( 



dq 



+ V(q) 



(62) 



where, in the last equality we have rescaled r — ► r//i. Eq. 62 is known as the 
Feynman-Kac Formula. 

Since the Partition Function is a trace over states, we must use boundary 
conditions such that the initial and final states are the same state, and to do 
a sum over all such states. In other terms, we must have periodic boundary 
conditions in imaginary time (PBC's),i.e., 



q{r) = q(r + (3) 



(63) 



Therefore a quantum mechanical system at finite temperature T can be de- 
scribed in terms of an equivalent system in classical statistical mechanics with 
Hamiltonian (or energy) 



(64) 



on a segment of length 1/T and obeying PBC's. This effectively means that the 
segment is actually a ring of length /3 = 1/T. 

Alternatively, upon inserting a complete set of eigenstates of the Hamilto- 
nian, it is easy to see that an arbitrary matrix element of the density matrix 
has the form 



(Q'\e- pH \q) =J2{q'\n){n\q)e-P E " 

oo 

= E e ~ /3£ "v>;(<z>«(<z) 



71=0 



(65) 



where {E n } are the eigenvalues of the Hamiltonian, E is the ground state 
energy and tpo(q) is the ground state wave function. 

Therefore, we can calculate both the ground state energy E and the ground 
state wave function from the density matrix and consequently from the (imagi- 
nary time) path integral. For example, from the identity 



Eq = — lim 



ilntre _/3ff 



(66) 
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we see that the ground state energy is given by 



If f fj 
E = — lim — In / Vq exp < — \ dr 

•P Jq(0)=q(P) ' I Jo 



(67) 



Mathematically, the imaginary time path integral is a better behaved object 
than its real time counterpart, since it is a sum of positive quantities, the sta- 
tistical weights. In contrast, the Feynman path integral (in real time) is a sum 
of phases and as such it is an ill-defined object. It is actually conditionally con- 
vergent and to make sense of it a number of convergence factors (or regulators) 
will have to be introduced. The effect of these convergence factors is actually an 
analytic continuation to imaginary time. We will encounter the same problem 
in the calculation of propagators. Thus, the imaginary time path integral, often 
referred to as the Euclidean path integral (as opposed to Minkowski), can be 
used to describe both a quantum system and a statistical mechanics system. 

Finally, we notice that at low temperatures T — ► 0, the Euclidean Path 
Integral can be approximated using methods similar to the ones we discussed 
for the (real time) Feynman Path Integral. The main difference is that we 
must sum over trajectories which are periodic in imaginary time with period 
(3 = 1/T. In practice this sum can only be done exactly for simple systems 
such as the harmonic oscillator, and for more general systems one has to resort 
to some form of perturbation theory. Here we will consider a physical system 
described by a dynamical variable q and a potential energy V(q) which has a 
minimum at q = 0. For simplicity we will take V(0) = and we will denote by 
tow 2 = V (0) (in other words, an effective harmonic oscillator). The partition 
function is given by the Euclidean path integral 



j Vq[l 



(68) 



where A E is the imaginary time, or Euclidean, version of the operator A, and 
it is given by 

Tfl d^ n 

A z = - ¥ J^ + V (g c (r)) (69) 

The functions this operator acts on obey periodic boundary conditions with 
period f3. Notice the important change in the sign of the term of the potential. 
Hence, once again we will need to compute a functional determinant, although 
the operator now acts on functions obeying periodic boundary conditions. In 
Physics 582 we will see that in the case of fermionic theories, the boundary 
conditions become antiperiodic. 

5.2.1 Computation of the Functional Determinant 

We will now do the computation of the determinant in i?( 2 ) . We will do the cal- 
culation in imaginary time and then we will carry out the analytic continuation 



12 



to real time. We want to compute 



D = Dct 



(70) 



subject to the requirement that the space of functions that the operator acts on 
obeys specific boundary conditions in (imaginary) time. We will be interested 
in two cases: (a) Vanishing Boundary Conditions (VBC's), which are useful 
to study quantum mechanics at T = 0, and (b) Periodic Boundary Conditions 
(PBC's) with period (3 = l/T. The approach is somewhat different in the two 
situations. 

1. Vanishing Boundary Conditions: 

This method is explained in detail in Sidney Coleman's book, Aspects of 
Symmetry. I will follow his approach closely. 

We define the (real) variable x = ^r. The range of x is the interval [0, L], 
with L = h(3/y/m. Let us consider the following eigenvalue problem for 
the Schrodinger operator — d 2 + W(x), i.e., 

(-d 2 + W(x)) ip(x) = X^(x) (71) 

subject to the boundary conditions V(0) = = 0. Formally, the 

determinant is given by 

D = ~[[\ n (72) 

n 

where {A„} is the spectrum of eigenvalues of the operator — d 2 + W(x) for 
a space of functions satisfying a given boundary condition. 

Let us define an auxiliary function ip\ (x) , with A a real number not neces- 
sarily in the spectrum of the operator, such that the following requirements 
are met: 

(a) tp\(x) is a solution of Eq. 71, and 

(b) tp\ obeys the initial conditions, tp\(0) = and d x ip\(0) = 1. 

It is easy to see that — d 2 + W(x) has an eigenvalue at A„ if and only 
if i/j\ n (L) = 0. (Because of this property this procedure is known as the 
Shooting Method.) Hence, the determinant D of Eq. 72 is equal to the 
product of the zeros of ip\(x) at x = L. 

Consider now two potentials and W^ 2 \ and the associated functions, 
ip^\x) and ip^\x). Let us show that 



Det 



-dP + WW(x)-\ 



£!m (73 ) 



The 1. h. s. of Eq. 73 is a meromorphic function of A in the complex plane, 
which has simple zeros at the eigenvalues of — d 2 + W^ x \x) and simple 
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poles at the eigenvalues of — d 2 + W^(x). Also, the 1. h. s. of Eq. 73 
approaches 1 as |A| — > oo, except along the positive real axis which is 
where the spectrum of eigenvalues of both operators is. Here we have 
assumed that the eigenvalues of the operators are non-degenerate, which 
is the general case. Similarly, the r. h. s. of Eq. 73 is also a meromorphic 
function of A, which has exactly the same zeros and the same poles as the 
1. h. s. . It also goes to 1 as |A| — > oo (again, except along the positive 
real axis), since the wave- functions ip\ are asymptotically plane waves in 
this limit . Therefore, the function formed by taking the ratio r. h. s. / 1. 
h. s. is an analytic function on the entire complex plane and it approaches 
1 as |A| — > oo. Then, general theorems of the Theory of Functions of a 
Complex Variable tell us that this function is equa7 to 1 everywhere. 

From these considerations we conclude that the following ratio is indepen- 
dent of W{x), 

Dct (-d 2 + W{x) - A) 



ML) 

We now define a constant Af such that 

Dct (-d 2 + W{xj) 



ML) 



Then, we can write 



Af [Dct (-d 2 + W)] 1/2 = [nhip {L)] 



-1/2 



(74) 



(75) 



(76) 



Thus we reduced the computation of the determinant (including the nor- 
malization constant) to finding the function ipo(L). For the case of the 
linear harmonic oscillator, this function is the solution of 

d 2 

'dx~ 2+m 



Mx) = o 



(77) 



with the initial conditions, Vo(0) = and tp' (0) = 1. The solution is 

1 



■ sinh(v / rowa;) 



(78) 



Hence, 



Z = M 



Det (--^ + mu; 2 ^ = [^o(£)P 1/2 (79) 



and we find 



irh 



sinh(/3w) 



-1/2 



(80) 



where we have used L == hfl/y/m. From this result we find that the 
ground state energy is 



Eq 



lim — — In Z 



~2~ 



(81) 
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as it should be. 

Finally, by means of an analytic continuation back to real time, we can 
use these results to find, for instance, the amplitude to return to the origin 
after some time T. Thus, for tf —U = T and qj = qi = 0, we get 



(0,T|0,0) = 
2. Periodic Boundary Conditions: 



sin(wT) 



-1/2 



(82) 



Periodic boundary conditions imply that the histories satisfy q(r) = q(r + 
/?). Hence, these functions can be expanded in a Fourier series of the form 



Q(r)= £ e^q n 



(83) 



where u> n = 2nn/(3. Since q(r) is real, we have the constraint q- n = q*. 
For such configurations (or histories) the action becomes 



dr 



m 

2h 2 \dr 



-V (0)q' 



2 

-^"(0)% 2 + fE[^ + ^"(°) 



2 



n>l 



The integration measure now is 
Vq[r 



y dq -q (IJrnq,, 



n>l 



2ir 



(84) 



(85) 



where Af is a normalization constant that will be discussed below. After 
doing the Gaussian integrals, the partition function becomes, 

1/2 



n 



i 



n 



i 



_n— — oo K 



(86) 

Formally, the infinite products that enter in this equation are divergent. 
The normalization constant N eliminates this divergence. This is an ex- 
ample of what is called a regularization. The regularized partition function 
is 

h 2 v"(o)~ 



z = 

Using the identity 



1 



h 2 (3 y/pv"(o) 

nfi+i 



n 



i + 



n>l 



sinha 



(87) 



(88) 
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we find 

Z = ; 

' ph I v"(o) 



z = 7 ; ( 89 ) 



2 sinh 



which is the standard partition function for a linear harmonic oscillator, 
see L. D. Landau and E. M. Lifshitz, Statistical Physics. 

5.3 Path Integrals for a Scalar Field Theory 

We will now develop the path-integral quantization picture for a scalar field 
theory. Our starting point will be the canonically quantized scalar field. As we 
saw before in canonical quantization the scalar field 4>(x) is an operator that acts 
on a Hilbert space of states. We will use the field representation, which is the 
analog of the conventional coordinate representation in Quantum Mechanics. 
Thus, the basis states are labelled by the field configuration at some fixed time 
xo, i.e., a set of states of the form { \{4>(x, xo)}} }. The field operator <f>[x, xq) 
acts trivially on these states, 

4>(x,x )\{(p(x,x )}) = <j>(x,x )\{<l>(x,xo)}} (90) 

The set of states { |{</>(x, %o)}) } is both complete and orthonormal. Complete- 
ness here means that these states span the entire Hilbert space. Consequently 
the identity operator T in the full Hilbert space can be expanded in a complete 
basis in the usual manner, which for this basis it means 

1 = J Vct>(x,x ) \{<t>(2,xo)})({<f>(2,xo)}\ (91) 

Notice that since the completeness condition involves a sum over all the states 
in the basis and since this basis is the set of field configurations at a given time 
xo, we will need to give a definition for integration measure which represents 
the sums over the field configurations. In this case there is a trivial definition, 

V<f>(2,x ) =Y[d<j)(x,x ) (92) 

x 

Likewise, othonormality of the basis states is the condition 

(<t>(x,x )\<t>'(x,xo)) = Y[S (<t>(x,x ) - 4>'{x,x Q )) (93) 

X 

Thus, we have a working definition of the Hilbert space for a real scalar field. 
Naturally, there are many other definitions of this Hilbert space and they are 
all equally good. 

We saw in the previous section that in canonical quantization the classical 
canonical momentum n(x, xo), defined as 

SjC- 

n(x, X ) = jjzi r = d (p(x, Xq) (94) 

d<9 <?Kx,xo) 



16 



becomes an operator that acts on the same Hilbcrt space as the field itself <j> 
does. The field (j> and the canonical momentum ir satisfy equal time canonical 
commutation relations (CCR) 

x,x ),tl(y,x ) =ih5 3 (x-y) (95) 

Here we will use the Lagrangian density for a real scalar field 

1 (8..^ 2 



£=7; (V) ~ V(4>) (96) 



It is a simple matter to generalize what follows below to more general cases, 
such as complex fields and/or several components. Let us also recall that the 
Hamiltonian for a scalar field is given by 



/Tft 2 1 / - -\ 2 
d 3 x — + -(y ( f,J +V(tf>) 



(97) 



For reasons that will become clear soon, it is convenient to add an extra term 
to the Lagrangian density of the scalar field, Eq. 96, of the form 

^source = J ( x ) 4>{x) (98) 

The field J(x) is called an external source and represents the effects of external 
sources on the scalar field. The field J(x) is the analog of external forces acting 
on a system of classical particles. Here we will always assume that the sources 
J(x) vanish both at spacial infinity (at all times) and everywhere in both the 
remote past and in the remote future, 

lim J(x, xq) = lim J(x, xq) = (99) 

|x|— >oo x ^±oo 

The total Lagrangian density is 

£(4>, J) = £ + 

^source 

(100) 

Notice that since the source J (a;) is generally a function of space and time, the 
Hamiltonian that follows from this Lagrangian is formally time-dependent. 

We will derive the path integral for this quantum field theory by following 
the same procedure we used for the case of a finite quantum mechanical system. 
Hence we begin by considering the amplitude 

j({<P(x,x )}\{cp'(y,y )})j (101) 

In other words, we want the amplitude in the background of the sources J{x). 
We will be interested in situations in which xq is in the remote future and yo is 
in the remote past. It turns out that this amplitude is intimately related to the 
computation of ground state (or vacuum) expectation values of time ordered 
products of field operators in the Heiscnberg representation 

G (N \x u ...,x N ) = (0\T[4>(x 1 )...4>(x N )}\0) (102) 
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which are the TV-point Green Functions. In particular the 2-point function 

G ( - 2 \x 1 ,...,x N ) = (0\T[4>(x 1 )4>(x 2 )]\0) (103) 

is known as the Feynman Propagator for this theory. We will see later on that 
all quantities of physical interest can be obtained from a suitable Green function 
of the type of Eq. 102. 

In Eq. 102 we have use the notation T[<f>{x\) . . . <^(xjv)] for the time-ordered 
product of Hcisenbcrg field operators. For any pair Hciscnbcrg of operators 
A(x,xq) and B(y,y ), (which commute for space-like separations) the time or- 
dered product is defined to be 

T[A(x, x )B(y, y )} = 0{x - y )A(x, x )B(y, y ) + 6(y - x )B(y, y )A(x, x Q ) 

(104) 

where 9{x) is the step (or Heaviside) function 

'(*)=(; if ;- ' d05) 

I otherwise 

This definition is generalized by induction to any number of operators. Notice 
that inside a time-ordered product the Heisenberg operators behave as numbers. 

Let us now recall the structure of the derivation that we gave of the path 
integral in Quantum Mechanics. We will paraphrase that derivation for this 
field theory. We considered an amplitude equivalent to Eq. 101, and realized 
that this amplitude is actually a matrix element of the evolution operator, 



dx' Q H(x' )/h 

r({<t>(x',x )}\{<j>'(y,yo)})j = ({<f>(x')}\Te J v« \W{v)}) (106) 

where T stands for the time ordering symbol (not temperature!), and H(x' ) is 
the time-dependent Hamiltonian: 

H(x ) = J d 3 x ^U 2 (x,x ) + i (V(f>(x,x )^ + V(<f>(x,x )) - J(x,x ) 

(107) 

We then partitioned the time interval in a large number of steps of width At 
and inserted a complete set of eigenstates of the field operator 0, since it plays 
the role of the coordinate. As it turned out, we also had to insert complete 
sets of eigenstates of the canonical momentum operator, which here means the 
operator Il(x). The result is the phase-space path integral 

, - / dx' 4>u - H{<t>, n) + j(t> 

A{mxo)}\{<i>'(y,yo)})j= V<pvile n Jy° 

Jb. c. 

(108) 

where b. c. indicates the boundary conditions required by the requirement that 
the initial and final states be \{(j)(x,xo)}} and \{<fr'{y,yo)}) respectively. 
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Exactly as in the case of the path integral for a particle, this theory has a 
Hamiltonian quadratic in the momenta H(x). Hence, we can further integrate 
out the field H(x), and obtain the Feynman path integral for the scalar field 
theory in the form of a sum over histories of field configurations: 



,({<j>(x,x )}\{<f>'(y,y )})j=Af f V<p eh 

Jb. c. 



-S(<j>,d^<j>,J) 



(109) 



where N is an (unimportant) normalization constant, and S{<j>, d^cj), J) is the 
action for a real scalar field (f>(x) coupled to a source J(x), 



S(<f>,d^,J) = J d 4 



I (d^f - V{4>) + H 



(110) 



5.4 Path Integrals and Green Functions 

In Quantum Field Theory we will be interested in calculating vacuum (i. e. 
ground state) expectation values of field operators at various space-time lo- 
cations. Thus, instead of the amplitude j({4>(x, xo)}\{(f>' (y,y )}} j we may be 
interested in a transition between an initial state, at yo — > — oo which is the 
vacuum state |0), i.e., the ground state of the scalar field in the absence of the 
source J (a;), and a final state at x — > oo which is also the vacuum state of the 
theory in the absence of sources. We will denote this matrix element by 

Z[J} = j(0\0)j (111) 

This matrix element is called the Vacuum Persistence Amplitude. 

Let us see now how the vacuum persistence amplitude is related to the 
Feynman path integral for a scalar field of Eq. 109. In order to do that we will 
assume that the source J(x) is "on" between times t < t' and that we watch 
the system on a much longer time interval T < t < t' < T' . For this interval, 
the amplitude is 

j({&(x,T')}\{$(x,T)})j = J Vcj>{x,t) Vcf>'(x,t') 

({$'(^r')}|{^(f,f)}>j<{^ , (^0}i>{^*)} J mx,t)}mx,T)}) 

(112) 

The matrix elements ({$'(£, T')}\{(j)'(x, t')}} and ({<j>(x, t)}\{$(x, T)}) are given 
by 

m 

({$'(x, t')}\w(x, f)}) * n [{&mmwm\ e-^ T ' 

n 

(113) 
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where we have introduced complete sets of eigenstates of the Hamilto- 

nian of the scalar field (without sources) and the corresponding wave functions, 
{*„($((*))}. 

We now analytically continue T along the positive imaginary time axis, and 
T' along the negative imaginary time axis, as shown in figure 3. After carrying 



Im t 








T 








Re t 




I V 






Figure 3: Analytic continuation 





out the analytic continuation, we find that the following identities (essentially, 
the Gell-Mann-Low Theorem) hold, 

lim e-^ T l h mx,t)}\{<$>{x,T)}) = *o[{</>}] %[{$>}} e~ iEot / h 

lim e iE ° T '/ H {{<f>\x,T')}\{t(x,t')}) = *o [{#'}] e lE ^' l h 

(114) 

and all other terms drop out in this limit provided the vacuum state |0) is 
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non-degenerate. Hence, in the same limit, we also find the following relation: 

lim lim ({n*,T) }mS , m 

T^+iooT'^-ioo exp[-iE (T' -T)/h] ^o[{<5>}} *o[{$'}] 

= J v<s>v<s>'%M(x,t')}} *o[{0(*,t)}M{^W)}IW2>*)}>j =j(o\o)j 

(115) 

Eq. 115 gives us a relation between the Feynman Path Integral and the vacuum 
persistence amplitude of the form 



Z[J} = j(0\0)j=M lim lim fv^e 



T' 

4 



d 4 x [c(<f>, d^) + J4>] 



(116) 

In other words, in this asymptotically long time limit, the amplitude of Eq. 101 
becomes identical to the vacuum persistence amplitude j(0|0)j, independently 
of the choice of the initial and final states. 

Hence we find a direct relation between the vacuum persistence function 
Z[J] and the Feynman Path Integral, given by Eq. 116. Notice that in this 
limit we can ignore the "hard" boundary condition and work instead with free 
boundary conditions. Or equivalently, physical properties become independent 
of the initial and final conditions placed. For these reasons, from now on we 
will write the simpler expression 

r % - [d 4 X [£(4>,d ll 4>) + J<t>] 
Z[J] = j (0\0) j =Af V<j)e h J (117) 

This is a very useful relation. We will see now that Z[J] is the generating 
function(al) of all the vacuum expectation values of time ordered products of 
fields, i.e. the Green functions. In particular, let us compute the expression 



1 S 2 Z[J] 



Z[0] 6J(x)6J(x') 



S 2 j(0\0) 



j 



2 



(0\T{<P(xMx')}\0) 



j=v <0|0) 5J(x)5J(x>) 

(118) 

Thus, the 2-point function, the Feynman propagator or Green function of the 
scalar field <fi(x), (0\T{<f>(x)<j>(x')] |0) becomes, 

(0\T[<f>(x)<f>(x')}\0) = ^Jv<f, <f>(x) 4>(x') cxp (^-S[cj>, (119) 

Similarly, the iV-point function (0|T[(/)(x 1 ) . . . <f)(x N )]\0) becomes 



(o\T[<f>( Xl )...<p(x N )}\o) = (-ih) N 1 f>N J<;<)|()> -' 



j=o 



(0|0) 5J(x{)...5J{x N ) 
= -^y J Vcj> (pixi) . . . <p{x N ) exp {^S[<j>, d^] 



(120) 
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where 



Z[0] = (0|0)= f V<j> expf^Sfad^] 



(121) 



Therefore, we find that the Path Integral always yields vacuum expectation 
values of time-ordered products of operators. The quantity Z[J] can thus be 
viewed as the generating functional of the Green functions of this theory. These 
are actually general results that hold for the path integrals of all theories. 

5.5 Path Integrals in Euclidean Space and Statistical Physics 

In the last section we saw how to relate the computation of transition amplitudes 
to path integrals in Minkowski space-time with specific boundary conditions 
dictated by the nature of the initial and final states. In particular we derived 
explicit expressions for the case of fixed boundary conditions. However we could 
have chosen other boundary conditions. For instance, for the amplitude to begin 
in any state at the initial time and to go back to the same state at the final 
time, but summing over all states. This is the same as to ask for the trace 



where PBC stands for periodic boundary conditions on some generally finite 
time interval t' — t, and T is the time-ordering symbol. 

Let us now carry the analytic continuation to imaginary time t — > —it, i.e. 
a Wick Rotation. Upon a Wick rotation the theory has Euclidean invariance, 
i.e., rotations and translations in D = d+ 1-dimensional space. Imaginary time 
plays the same role as the other d spacial dimensions. Hereafter we will denote 
imaginary time by xd, and all vectors will have indices \i that run from 1 to D. 

We will consider two cases: infinite imaginary time interval, and finite imag- 
inary time interval. 

5.5.1 Infinite Imaginary Time Interval 

In this case the path integral becomes 



where D is the total number of space-time dimensions. Here we are discussing 
the case D = 4, but it is obviously valid more generally. Here He is the Eu- 
clidean Lagrangian 




(122) 




(123) 



Ce 



^(do0) 2 + i(V0) 2 + ^) 



(124) 
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Figure 4: Periodic boundary conditions wraps space-time into a cylinder. 

The path integral of Eq. 123 has two interpretations. One is simply the infinite 
time limit (in imaginary time) and therefore it must be identical to the vac- 
uum persistence amplitude j(0|0) J . The only difference is that from here we 
get all the TV-point functions in Euclidean space-time (i.e., imaginary time). 
Therefore, the relativistic interval is 

t 2 - x 2 -t 2 - x 2 < (125) 

which is always space-like. Hence, with this procedure we will get the Green 
functions for space-like separations of its arguments. To get to time-like sepa- 
rations we will need to do an analytic continuation. This we will do later on. 
The second interpretation is that the path integral of Eq. 123 is the partition 
function of a system in Classical Statistical Mechanics in D dimensions with 
energy density (divided by T) equal to Ce — J<f>- This will turn out to be a very 
useful connection (both ways!). 

5.5.2 Finite Imaginary Time Interval 

In this case we have 

< x = t < 13 = 1/T (126) 

where T will be interpreted as the temperature. Indeed, in this case the path 
integral is 

Z'[0]=1ie~^ H (127) 

and we are effectively looking at a problem of the same Quantum Field Theory 
but at finite temperature T — 1/(3. The path integral is once again the partition 
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function but of a system in Quantum Statistical Physics! The partition function 
thus is (I set h~ = 1) 



r dT(jC E -J<f>) 

/ T>4> e Jo 



Z'[J] = J V(f> e Jo (128) 

where the field <j>(x, t) obeys periodic boundary conditions in imaginary time, 

(j>(x,T)=(f>(x,T + f3) (129) 

This boundary condition will hold for all bosonic theories. Theories with fermions 
obey instead, as we will see later on, anti-periodic boundary conditions. 

Hence, Quantum Field Theory at finite temperature T is just Quantum 
Field Theory on an Euclidean space-time which is periodic (and finite) in one 
direction, imaginary time. In other words, we have wrapped ( "compactified" ) 
Euclidean space-time into a cylinder with perimeter (circumference) (3 = 1/T 
(in units of h = ks = 1). 

Clearly, the Green functions in imaginary time (which we will call the Eu- 
clidean Correlation functions) are given by 

XN yl \ Tl 

" " ! = (4>{x 1 )...</,{x N )) (130) 



Z'[J) 5J{xi)...J{x N ) 



which are just the correlation functions in the equivalent problem in Statistical 
Mechanics. Upon analytic continuation the Euclidean Correlation Functions 
(4>(xi) . . . 4>{xn)) and the iV-point functions of the QFT are related by 

(<P( Xl ) . . . cj>{x N )) <- (ih) N (Q\T<P( Xl ) . . . 4>{x N )\Q) (131) 



For the case of a QFT at finite temperature T, the path integral yields the 
correlation functions of the Hciscnberg field operators in imaginary time. These 
correlation functions are often called the thermal Green functions. They are 
functions of the spatial positions of the fields, x\, . . . ,xn and of their imaginary 
time coordinates, Xdi, ■ ■ ■ , xdn (here Xd = t). To obtain the correlation func- 
tions as a function of the real time coordinates x i, ■ ■ ■ , x 0N at finite temperature 
T it is necessary to do an analytic continuation. We will discuss how this is done 
later on. 



5.6 Path Integrals for the Free Scalar Field 

We will consider now the case of a Free Scalar Field. We will carry our discussion 
in Euclidean Space-Time (i.e., in imaginary time), and we will do the relevant 
analytic continuation back to real time at the end of the calculation. 

The Euclidean Lagrangian Ce for a free field 4> coupled to a source J is 

Ce = \ (V^) 2 + \™ 2 tf - H (132) 
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where we are using the notation 

(V^) 2 = V^W<A (133) 

Here the index is /j, — 1, . . . ,D for an Euclidean space-time of D = d + 1 
dimensions. For the most part we will be interested in the case of d = 3 and 
Euclidean space has four dimensions. Notice the way the Euclidean space-time 
indices are placed in Eq. 133. This is not a misprint! 

We will compute the Euclidean Path Integral (or Partition Function) Ze[J] 
exactly. The Euclidean Path Integral for a free field has the form 



Z E [J]=M Jvcj) e 



\ (v^) 2 + \™ 2 4> 2 - H 



(134) 



In Classical Statistical Mechanics this theory is known as the Gaussian Model. 

In what follows I will assume that the boundary conditions of the field (j> 
(and the source J) at infinity are either vanishing or periodic, and that the 
source J also either vanishes at spatial infinity or is periodic. With these as- 
sumptions all terms which are total derivatives drop out identically Therefore, 
upon an integration by parts and after dropping boundary terms, the Euclidean 
Lagrangian becomes 

C E = \<t> [- XJ 2 +m 2 ] <}>- 3<\> (135) 

This path integral can be calculated exactly because the action is a quadratic 
form of the field <fi. It has terms which are quadratic (or, rather bilinear) in <f> 
and a term linear in <j>, the source term. By means of the following shift of the 
field 

tj>{x) = 4>{x) + £(x) (136) 

the Lagrangian becomes 
Ce = \i> [- V 2 +m 2 ] <j)-Jc^ 

(137) 



Hence, we can decouple the source J(x) by requiring that the shift (f> be such 
that the terms linear in £ cancel each other exactly. This requirement leads to 
the condition that (j> be the solution of the following partial differential equation 

[- v 2 +m 2 ] 4> = J{x) (138) 

Equivalently we can write <j) is terms of J(x) through the action of the inverse 
of the operator — y 2 +m 2 , 

4> - -4- 2 J (139) 
— v + m 



25 



The solution of Eq. 138 is 



<p(x) = J d D x' G E (x - x 1 ) J(x') (140) 

where 

G*(x-x') = {x\ _J +m2 \x>) (141) 

is the Green function of the linear partial differential operator — v 2 +m 2 , i.e., 
Gq(x — x') is the solution of 

[- V l +m 2 ] G £ (x - x') = S D (x- x') (142) 

In terms of Gq(x — x'), the terms of the shifted action become, 

J d D x (±4>{x) [- v 2 +m 2 ] fa) - Jfa)) = J d D x (-±${x)J{x)) 
= - 1 - Jd D x d D x' J(x) G%(x - x') J(x') 

(143) 

Therefore the path integral Ze[J], defined in Eq. 134, is given by 

\ ( d D xd D x' J{x)GZ{x-x') J(x') 
Ze[J] = Z e [0] e 2 J (144) 

where Ze[0] is 

-\ f A^)[- V 2 W]^) 

V£, e 2 J (145) 



,-[0]=/ 



Eq. 144 shows that, after the decoupling, Ze[J] is a product of two factors: 
(a) a factor that is function of a bilinear form in the source J, and (b) a path 
integral, ^_e[0], that is independent of the sources. 

5.6.1 Calculation of Z E [0] 

The path integral Ze[0] is analogous to the fluctuation factor that we found in 
the path integral for a harmonic oscillator in elementary quantum mechanics. 
There we saw that the analogous factor could be written as a determinant of a 
differential operator, the kernel of the bilinear form that entered in the action. 
The same result holds here as well. The only difference is that the kernel is now 
the partial differential operator A = — \/ 2 +m 2 whereas in Quantum Mechanics 
is an ordinary differential operator. Still, the operator A has a set of eigenstates 
{^n(x)} which, once the boundary conditions in space-time are specified, are 
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both complete and orthonormal, and the associated spectrum of eigenvalues A n 
is 

[- v 2 +™ 2 ] *„(&) = A n $ n (x) 
J dPx tf n (ar) * m (x) = S n , m 
2 tt„(aO *„(x') = - x') 

n 

(146) 

Hence, once again we can expand the field <f>(x) in the complete set {^ n (x)} , 

<P(x) =^c„ *„(x) (147) 

n 

The set of field configurations is thus parametrized by the coefficients {c„}. 
The action now becomes, 

S= I ' d D x C E {4>, dt) = \Y, AnC " ( 148 ) 

n 

Thus, up to a normalization factor, we find that 2"_e[0] is given by 

Ze[0] =U A n 1/2 = ( Dct [- V 2 W])- 1/2 (149) 

n 

Once again, we have reduced the calculation of Ze [0] to the computation of the 
determinant of a differential operator, Det [— V 2 +m 2 ~\ . 

In chapter 8 we will discuss efficient methods to compute such determinants. 
For the moment it will be sufficient to notice that there is a simple, but formal, 
way to compute this determinant. First, we notice that if we are interested in 
the behavior of of an infinite system at T = 0, the cigenstates of the operator 
— v 2 +m 2 are simply suitably normalized plane waves. Let L be the linear size 
of the system (L — > oo). Then the eigenfunctions are labeled by a D-dimensional 
momentum (with fi = 0, 1, . . . , d) 

^ = j2^ eiP " X " (150) 

with eigenvalues, 

A p =p 2 +m 2 (151) 
Hence the logarithm of determinant is 

In Det [- v 2 +m 2 ] = Tr In [- \j 2 +m 2 ] 

= ^2 m (p 2 + m2 ) 

p 

(152) 
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where V = L D is the volume of Euclidean space-time. Hence, (the logarithm 
of) Z E [0] is 

InZ B [0] = -\ J Hp 2 + m 2 ) (153) 

This expression is has two singularities: 

1. It diverges as V — > oo. This infrared (IR) singularity actually is not a 
problem since lnZefO] should be extensive. In other words, this is how it 
should behave. 



2. The integral diverges at large momenta unless there is an upper bound (or 
cutoff) for the allowed momenta. This is an ultraviolet (UV) singularity. 
It has the same origin of the UV divergence of the ground state energy. 

In fact Ze[0] is closely related to the ground state (vacuum) energy since 

Z E [0] = lim V e-P En ~ e~^ E ° + ... (154) 

(3— >oo ^ — 4 

n 

Thus, 

Eo = - lim I In Z E [0] = \h d J ^ ln(p 2 + m 2 ) (155) 

where L d is the volume of space, i. e. V = L d (3. Notice that Eq. 155 is UV 
divergent. Below in this chapter we will discuss how to compute expressions of 
the form of Eq. 155. 



5.6.2 Green Functions 



A number of interesting results are found immediately by direct inspection of Eq. 
144. We can easily see that, once we set J = 0, the Green function G^(x — x') 



G^(x-x') = (x\ 



1 



\x>) 



1 — v 2 +m 2 

is equal to the 2-point correlation function for this theory (at J = 0) 

1 S 2 Z E [J] 



(<t>{x)<t>(x')) = 



Z E [0] SJ{x)6J{x') 



= G° E (x-x') 



(156) 



(157) 



j=o 



Likewise we find that, for a free field theory, the TV-point function (<p(xi) . . . 4>(xn)) 
is equal to 



(<j){xi)...<t>(x N )) 



1 



S N Z E [J] 



Z E [0] 8J{xi)...8J{x N ) 



= (<j){x\)4>{x2)) ■ ■ ■ {(I)(xn-i)<P(xn)) + permutations 



(158) 
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Therefore for a free field, up to permutations of the coordinates x\, . . . ,Xn, 
the iV-point functions factorize into products of 2-point functions. Hence, N 
must be an positive even integer. This result, Eq. 158, which we derived in the 
context of a theory for a free scalar field, is actually much more general. It is 
known as Wick's Theorem. It applies to all free theories i.e., theories whose 
Lagrangians are bilinear in the fields, it is independent of the statistics and on 
whether there is relativistic invariance or not. The only caveat is that, as we 
will see later on, for the case of fermionic theories there is a sign associated with 
each term of this sum. It is easy to see that, for N = 2k, the total number of 
terms in the sum is 

(2fc)! 
2 k k\ 

Each factor of a 2-point function (<fi(xi)(j)(x2)), i.e. a free propagator, and it is 
also called a contraction. It also common to use the notation 



(2fc-l)(2fc-3) = (159) 



(0(xi)0(x 2 )) = <f>( Xl )<l>(x2) (160) 
to denote a contraction or propagator. 

5.6.3 Calculation of the Propagator 

We will now calculate the 2-point function, or propagator, Gq(x— x') for infinite 
Euclidean space. This is the case of interest in QFT at T = 0. Later on we will 
do the calculation of the propagator at finite temperature. 

Eq. 142 tells us that Gq{x — x') is the Green function of the operator — y 2 
+TO 2 . We will use Fourier transform methods and write Gq(x — x') in the form 

G*(x --') = / (0 Gf?(p) O (161) 

which is a solution of Eq. 142 if 

<#<■» = ^ (""J 

Therefore the Green function in real (Euclidean!) space is the integral 

G$(x -x')= / 44n 5 ^ (163) 

We will often encounter integrals of this type and for that reason we will do this 
one in some detail. We begin by using the identity 

A 

1 1 f°° 

A = 2 J, daG (164) 

where A > is a positive real number. The variable a is called a Feynman- 
Schwinger parameter. 
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We now choose A = p 2 +m 2 , and substitute this expression back in Eq. 163, 
which takes the form 



, 1 f°° f d D p -^(p 2 + m 2 )+ip ll {x^-x^) 



G*(x-x') = ^ da i-^e 2- — (165) 



The momentum integrand is now a Gaussian, and the integral can be calculated 
by a shift of the integration variables p^, i.e., by completing squares 

f (p 2 + m 2 ) - ip,(x, - = \ (V5p„ - 2 - \ 2 (166) 

and by using the Gaussian integral 



J (2n) D 



2 

D/2 



(2ira)- u/2 (167) 



(27T) 

After all of this is done, we find the formula 

\X-X'\ 2 1 2 
1 /-OO —J ! 777 (V 

w^L daa e 20 2 

(168) 

Let us now define a rescaling of the variable a, 

a = Xt (169) 

by which 

lX ~ ^ + ^m 2 a = ^ ~ ^ + l -m 2 \t (170) 



2a 2 2At 

We choose 



\x — x n 



A = J L (171) 

m 



With this choice, the exponent becomes 

\x — x'\ 2 1 9 m\x — x'\ ( 1 



2a + 2 ma =^r^[ t+ -t) ^ 
After this final change of variables, we find that the Green function is 

G °^"') = (^72 (|^7|j 2 ^-.(^-^'D (173) 

where K u (z) is the Bessel function (of imaginary argument, see e.g. Gradshteyn 
and Ryzhik) which has the integral representation 

z ( 1 



i c°° — — i 6 i I 

o / dtt v ~ 1 e 2 V V (174) 
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Figure 5: Behaviors of the Euclidean Green Function. 

with v = y — 1, and z = m|a; — 

There are two interesting regimes: (a) long distances, m\x — x'\ ^> 1 , and 
(b) short distances, m|a; — x'| <c 1. 

1. Long Distance Behavior. 

In this regime, z = m|x — x'\ 3> 1, a saddle-point calculation shows that 
the Bessel Function K v [z) has the asymptotic behavior, 

K v {z)K^e- z [1 + 0(1/*)] (175) 
Thus, in this regime the Euclidean Green function behaves like 
^j2m D - 2 e -m\x-x'\ 



G${x-x>) 



D - 1 



1 + H rr) 

mix — x'\ 



(2tt) d / 2 (m\x-x'\) 2 

(176) 



Therefore, at long distances, the Euclidean (or imaginary time) Green 
function has an exponential decay with distance (and imaginary time). 
The length scale for this decay is 1/m, which is natural since it is the only 
quantity with units of length in the theory. In real time, and in conven- 
tional units, this length scale is just the Compton wavelength, h/mc. In 
Statistical Physics this length scale is known as the correlation length £. 

2. Short Distance Behavior: 

In this regime we must use the behavior of the Bessel function for small 
values of the argument, 

KAz)K {jTf + 0{l/z- 2 ) (177) 
The Green function now behaves instead like, 

where . . . are terms that vanish as m\x — x'\ — > 0. Notice that the leading 
term is independent of the mass m. This is the behavior of the massless 
theory. 
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5.6.4 Behavior of the Green function in Minkowski space 



We must now address the issue of the behavior of the Green function in real 
time. This means that now we must do the analytic continuation back to real 
time x . In what follows we will set D = 4. 

Let us recall that in going from Minkowski to Euclidean space we contin- 
ued xo — ► — 1x4. There is also factor of i difference in the definition of the 
Green function. Thus, the Green function in Minkowski space Gq{x — x') is the 
expression that results from 

Go(x-x') = iG*(x-x')\ x ^. xo (179) 

The relativistic interval s is given by 

s 2 = {t-t'f - (x-x'f (180) 

The Euclidean interval (length) \x — x'\, and the relativistic interval s are related 

by 

\X - X'\ = y/(x - X 1 ) 2 -> (181) 

Therefore, in D = 4 space-time dimensions, the Minkowski space propagator is 

7 77-1 / 

G (x-x') = -l^-^K x {m^^) (182) 
we will need the asymptotic behavior of the Bessel function Ki(z), 



Ki(z) = 
K x (z) = 




, for z > 1 
, for z < 1 



(183) 

where C — 0.577215 ... is the Euler constant. Let us examine the behavior 
of Eq. 182 in the regimes: (a) space-like, s 2 < 0, and (b) time-like, s 2 > 0, 
intervals. 

1. {x-x') 2 =s 2 <0 

This is the space- like domain. By inspecting Eq. 182 we see that for space- 
like separations, the factor ^/—s 2 is a positive real number. Consequently 
the argument of the Bessel function is real (and positive), and the Green 
function is pure imaginary. In particular we see that, for s 2 < the 
Minkowski Green function is essentially the Euclidean Green function, 

Gq(x-x') = iG^{x-x') , fors 2 <0 (184) 
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Figure 6: Behaviors of the Green Function in Minkowski space-time. 



Hence, for s 2 < we have the behaviors, 



Gq(x — x 1 ) rts i- 



m 



-mi 



4tt 2 



(mV-s 2 ) 



3/2 



Go(x-x') 



47T 2 (-S 2 ) 



for m\ 
for m\ 



2. (x - x') 2 = s 2 > 

This is the time-like domain. The analytic continuation yields 



Go(x-x') 



47T 2 V / S 2 " 



Ki (imvr) 



s 2 » 1 

~^ < 1 

(185) 

(186) 



For pure imaginary arguments, the Bessel function K 1 (iz) is the analytic 
continuation of the Hankel function, K\(iz) — — f i?} 1 (— z) . This func- 
tion is oscillatory for large values of its argument. Indeed, we now get the 
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behaviors 



G (x-x') 
G (x - x') 



4tt 2 



3/2 



, for mv? 3> 1 



, for inVs^ <C 1 



(187) 



Notice that, up to a factor of i, the short distance behavior is the same 
for both time-like and space like separations. The main difference is that 
at large time-like separations we get an oscillatory behavior instead of an 
exponential decay. The length scale of the oscillations is, once again, set 
by the only scale in the theory, the Compton wavelength. 



5.6.5 Exponential decays and mass gaps 

The exponential decay at long space-like separations (and the oscillatory behav- 
ior at long time-like separations) is not a peculiarity of the free field theory. It 
is a general consequence of the existence of a mass gap in the spectrum. We can 
see that by considering the 2-point function of a generic theory, for simplicity 
in imaginary time. The 2-point function is 

G (2 \x-x" ,t-t') = (0\Tj>(x,T)4>(x" ,t')|0) (188) 

where T is the imaginary time-ordering operator. 

The Heisenberg representation of the operator <f> in imaginary time is (h = 1) 

0{x, t) = e Hr fa, 0) e~ HT (189) 

Hence, we can write the 2-point function as 

G^{x-x',t-t') = 
= 0(t - T')(0\e HT 0{x,O) e -#( T - T ') 4>(x',0) e- HT '\0) 
+ 0(t' -T)(0\e HT ' 0(f ,0) e -#( T '- r ) 4>(x 1 0)e- HT \0) 
= 6{t - t') eM T - T ') (O|0(x,O) e -#( r - r ') 0)|0) 
+ 9(t' - t) e E ^ T ' ~ T ) (O|0(x',O) e' 11 ^' ~ T "> #e,0)|0) 

(190) 

We now insert a complete set of eigenstates {|n)} of the Hamiltonian H, 



with eigenvalues {E n }. The 2-point function now reads, 

G (2) (x-x' ,t-t') = 
= 9(t-t') (0\4>{x,0)\n) (n|0(x',O)|O) e ~{ E n ~ Eq){t - t') 

n 

+ 6(t'-t) W(£',0)|n) (n|0(x,O)|O) e ~( E n ~ E )(t' - t) 



Since 



and that 



4>(x,0) = e iP -*4>(0,0)e- iP -* 



(191) 

(192) 
(193) 



P\0) = 0, P\n) = P n , 
where P n is the linear momentum of state |n), we can write 

(0|<Hx,0)|n) (n|^(f',0)|0) = |(O|0(O, 0)|n)| 2 e~ iPn ' f ') (194) 
using the above expressions we can write the expressions in Eq. 191 in the form 
G (2 \x-x" ,r-r') = £ IW(0,0)|n)| 2 

n 

+6>(r' - T)e~ iPn ' ( f ' - ^e"^" ~ So )( T ' ~ r ) 

(195) 

Thus, at equal positions, x = x", we obtain the simpler expression in the imag- 
inary time interval r — t' 

GW(0,t-t') = ^|(0|^(f,0)|n)| 2 x e -(-En-^o)|r-r'| 

n 

(196) 

In the limit of large imaginary time separation, |t — t'| — > oo, there is always 
a largest non-vanishing term in the sums. This is the term for the state |no) 
that mixes with the vacuum state |0) through the field operator <j>, and with the 
lowest excitation energy or mass gap E na — Eq. Hence, for |r — r'| — > oo, the 
2-point function decays exponentially like, 



G (2) (0,r-r') ~ |(O|0(f,O)|n o )| 2 x e -(^™o - ^o)|r - r'| 



(197) 



35 



Therefore, if the spectrum has a gap, the correlation functions (or Green func- 
tions) decay exponentially in imaginary time. In real time we get an oscillatory 
behavior. This is a very general result. Finally, notice that Lorentz invariance 
in Minkowski space-time (real time) implies rotational (Euclidean) invariance 
in imaginary time. Hence, exponential decay in imaginary times, at equal posi- 
tions, must imply (in general) exponential decay in real space at equal imaginary 
times (an laos in real time since the time difference in this case vanishes). Thus, 
in a Lorentz invariant system the propagator at space-like separations is always 
equal to the propagator in imaginary time. 
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